The "Majoranon" and how to reahze it in a tabletop experiment 



O 

(N 

o 

O: 

oo 



> 

o 
o 

d 

(N 



X: 



Changsuk Noh,^'EI B. M. Rodri'gucz-Lara,^'^ and Dimitris G. Angelakis^' -'^■[ll 

^Centre for Quantum Technologies, National University of Singapore, 2 Science Drive 3, Singapore 117542. 
^Instituto Nacional de Astrofisica, Optica y Electronica 
Calle Luis Enrique Erro No. 1, Sta. Ma. Tonantzintla, Pue. CP 72840, Mexico 
Science Department, Technical University of Crete, Chania, Crete, Greece, 73100 

We introduce the term Majoranon to describe particles that obey the Majorana equation, which 
are difTerent from tire Majorana fermions widely studied in various physical systems. A general 
procedure to simulate the corresponding Majoranon dynamics, based on a decomposition of the 
Majorana equation into two Dirac equations, is described in detail. It allows the simulation of the 
two-component chiral spinors, the building blocks of modern gauge theories, in the laboratory with 
current technology. Specifically, a Majoranon in one spatial dimension can be simulated with a 
single qubit plus a continuous degree of freedom, for example a single trapped ion. Interestingly, the 
dynamics of a Majoranon deviates most clearly from that of a Dirac particle in the rest frame, in 
which the continuous variable is redundant, making a possible laboratory implementation feasible 
with existing set ups. 
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Introduction - Quantum simulation was originally 
envisioned as an approach to study complex quantum 
systems that are difficult to understand using conven- 
tional methods [J However, it has recently been 
realized that the concept can also be used to engineer 
quantum dynamics not readily accessible in naturally 
occurring physical systems; e.g., elementary particles in 
lower spacetime dimensions. A notable example is an ex- 
perimental demonstration of quantum simulation of the 
Dirac equation [3, Q in 1-1-1 dimensional spacetime using 
trapped ions Q. There have been other proposals and a 
few experimental demonstrations regarding the quantum 
simulation of rclativistic equations and phenomena @- 
including a recent proposal to simulate the Majorana 
equation [0]. 

The Majorana equation is a Lorentz covariant equation 
discovered by Majorana as a generalization of the Dirac 
equation [13 Ell: 



(1) 



where we have set h = c = 1 for convenience. The gamma 
matrices 7'^ obey {7^,7''} = ^g^^"; the symbol V'c stands 
for charge conjugation of the spinor i/;: ij^c = (7(7°)-^ 1/;*, 
where C obeys C{YYC-^ = -7'^, = C'^ and 
C"^ = — C; g'^" is the standard metric in the vacuum [l^ . 
To avoid confusion, let us discuss here the nomenclature 
found in the literature. A field obeying the Majorana 
condition, i/jc = "01 is called a Majorana field (fermion, 
particle) , in which case Eq. ([T|) reduces to the Dirac equa- 
tion. The equation itself is sometimes called the Majo- 
rana equation when the field obeys the condition [l3j , but 
we will reserve the term Majorana equation for Eq. ([T]) 
when the condition does not hold and call the hypothet- 
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ical particles that obey it Majoranons. We would like to 
emphasize that Majoranons are thus different to Majo- 
rana fermions that are widely studied in the literature. 

The Majorana equation has received a renewed inter- 
est since the discovery of finite masses of neutrinos. A 
neutrino is a neutral fermion; therefore it could be either 
a Majorana fermion, in which case it is its own antipar- 
ticle, or a Dirac particle, in which case a neutrino will 
be different to an anti-neutrino. For this reason, most 
studies on Eq. ([T]) are limited to Majorana fermions and 
the Majorana equation by itself seems to be of academic 
interest only. However, Casanova et al. have recently 
proposed a general procedure to implement non-physical 
operations such as charge conjugation and time reversal, 
allowing for an experimental study of the equation Q. 

Closely related to the Majorana equation are the so- 
called two-component Majorana equations (TCMEs). In 
1-1-1 and 2-1-1 D, it turns out that one of the two TCMEs 
is equivalent to the Majorana equation, whereas in 3-1-1 
D both of them merge into the Majorana equation. The 
importance of the TCMEs is that they are obeyed by the 
building blocks of the modern gauge theory, chiral spinors 
[3, [1^ . In this work, we show that it is possible to de- 
compose one of the equations, say for a right-chiral field, 
into two Hamiltonian equations that can be simulated. 
Our method differs from the earlier proposal 0], which 
depends on enlarging the Hilbert space of the system to 
turn the Majorana equation into a higher-dimensional 
Dirac equation. An advantage is that the size of the 
Hilbert space required is smaller, but a tradeoff exists, 
in that a complete reconstruction of the spinor, includ- 
ing the continuous mode, is required. These points are 
discussed further below. 

This paper is organized as follows. We start by giving 
a short summary on quantum simulation of the Majo- 
rana equation based on Hilbert space expansion men- 
tioned above, followed by a brief introduction to TCMEs 
for non-specialists. We show how the equations can be 
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motivated by trying to introduce Lorentz invariant mass 
terms in the Weyl equations and describe how rel- 
ativistic equations for higher-dimensional spinors can be 
constructed from these equations. We then prove that 
a general fermion obeying what we will call a Dirac- 
Majorana equation can be decomposed into two Majo- 
rana particles [l^, [3 [13 j which implies that such an 
equation can be physically realized. Wc end by describ- 
ing a general simulation procedure based on this decom- 
position. 

Quantum simulation by Hilbert space expan- 
sion - As shown in it is possible to implement un- 
physical operations such as charge conjugation, complex 
conjugation, and time reversal, via mapping to an en- 
larged Hilbert space. It is thus possible to simulate the 
Majorana equation which in 1+1 D reads 

idtip = a^Pxi' - imayil:* , (2) 

where V' = € C so V' G C^- To see how, 

first note that this equation can be rewritten as 

idt {ip + ip*) = (TxPx (-0 + "0*) + (0 - , (3) 
idt {ip - i>*) = (JxPx (0 - 0*) - imay (-0 + 0*) ■ (4) 

Now, by mapping to an extended Hilbert space such 
that * = (Re(V'(^)),Re(0(2)),Ini(0(i)),Im(0(2))), ^ = 
Af^f, M = (l2,?l2), the Majorana equation can be writ- 
ten as the Schrodinger equation, idt'^ = Hm^, where 



Hm = (I2 ® a-x)px - m{ax ® a-y)- 



(5) 



This Hamiltonian, describing the Dirac equation, can be 
implemented with two trapped ions. 

Two-component formalism - Here we give a brief 
introduction to the two component formalism of chiral 
spinors and show how they are related to the usual four- 
component Dirac equation and the Majorana equation 
[3 H, [13 • Being the smallest irreducible representa- 
tions of the Lorentz group, the chiral (or Weyl) spinors 
are the basic building blocks of modern gauge theories 
[1^ . Moreover, chirality plays an important role in the 
Standard Model due to the chiral nature of electroweak 
interaction, i.e. the electroweak interaction distinguishes 
between left- and right-chiral fields. 

Let us start with the description of masslcss fcrmions 
by the Weyl equations: 



(dt-a- V)0L = 0, 
[dt+a- V)0fl, = 0, 



(6) 
(7) 



where a ~ {a^, cTy^az) arc the Pauli matrices. The Weyl 
spinors ipji and ip^ transform according to the two dimen- 
sional irreducible representation of the Lorentz group and 
its complex conjugate, respectively, and are called two- 
component chiral spinors. In the massless case, the chi- 
rality coincides with the helicity defined as the projection 
of a particle's spin along its momentum. 

The Icft-chiral Weyl spinor was believed to describe the 
masslcss neutrinos and had been incorporated in the old 



Standard Model before the discovery of non-zero neutrino 
masses. To introduce mass terms in the Weyl equations, 
one must be careful about the Lorentz covariance of the 
equation. In order to keep the equations covariant, one 
should note that me~^0*, where is either 0^ or ipn 
and e = iay, transforms like the differential part of the 
corresponding Weyl equation [l^ . Therefore it is possi- 
ble to write the Lorentz invariant equations for massive 
fcrmions as 

i{dt-a-V)i'L{x)-me^jjl{x)=0, (8a) 
iidt + a ■ V)i'R{x) + meV4(a;) = 0, (8b) 

called the Icft-chiral and right-chiral two-component Ma- 
jorana equations. These spinors form the building blocks 
of higher-dimensional spinors such as the Dirac or Ma- 
jorana fields. They do not correspond to any massive 
particles by themselves, however, and a direct observa- 
tion of their dynamics is ordinarily impossible. In fact, 
Eqs. (|H1) is equivalent to the Majorana equation, which 
can be easily seen by constructing the four-component 
spinor (0_r,0l)- In 1+1 or 2-1-1 D, the situation simpli- 
fies because a Majoranon can be represented by a two- 
component spinor and it can be shown, by explicit con- 
struction of 7'' in terms of Pauli matrices, that one of 
the two TCMEs is equivalent to the Majorana equation. 

It is also possible to build a Majorana field from either 
of the chiral fields. To see this, for 0/, say, note that c^'l 
obeys the equation 

i{dt + a ■V)<i%Pl{x) + meeipL{x)=Q. (9) 

That is, e?/'2 behaves like a right-chiral two-component 
spinor. Thus noticing that Eq. (|8^) mixes the left- and 
right-chiral spinors, let us define a four-component spinor 



which obeys the Dirac equation, 

il^dt,-^M - m-^M = 0, 



(10) 



where are the Dirac matrices in the chiral representa- 
tion. 

Wc have thus shown that a four-component spinor 
obeying the Dirac equation can be constructed from a 
two-component Icft-chiral spinor. Using the properties 
of the Dirac matrices it is possible to show that 'i'M has 
a special property 



(11) 



where 



c = i^'^y 

— zcr. 



in the chiral representation. The subscript c denotes 
charge conjugation and the (Majorana) condition states 
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that the spinor is invariant under charge conjugation. A 
Majorana particle is therefore necessarily charge neutral, 
which is why such a particle can be its own antiparticle. 

To build up a Dirac field from two-component spinors 
we need two independent left-chiral spinors, say ipL and 
XL- Then {eip1,XL)'^ obeys the Dirac equation, but not 
the Majorana condition. 

The fact that the TCME can be converted into the 
four-component Dirac equation suggests that a two- 
component Majoranon is equivalent to a four-component 
Majorana fermion. The resulting equation can then be 
simulated with a single trapped ion using the method 
proposed in Furthermore, if we restrict the spatial 
dimension to 1 it is possible to write down the equa- 
tion of motion in the form of interacting qubits coupled 
to a phonon mode, where only one qubit is coupled to 
the phonons. This fact has already been noticed in Q 
through different argument and allows quantum simula- 
tion of the Majorana equation using trapped two-level 
ions. Our argument here can thus be seen as an alter- 
native derivation of the result, starting from the two- 
component formalism. The salient point is that in 1+1 
D a Majoranon is described by a two-component spinor 
obeying either of Eqs. ^ with 1 spatial degree of free- 
dom and therefore it should be possible to convert this 
to a Hamiltonian dynamics by constructing a Majorana 
fermion as shown earlier. To see this explicitly, con- 
sider the Dirac equation for a Majoranon confined to 
move in the z-direction. In the Majorana representa- 
tion = o-y ® (Ta;, = iay (g) ay and 5'^^ = ^E'M, and 
Eq. becomes 

idt-^M = 1 ® (7,iidz)'i'M + m{ay ® a^)^M- (12) 

The Hamiltonian for this Schrodinger equation is equiva- 
lent to ([5]) up to an unitary transformation and thus can 
be readily simulated with trapped ions as shown explic- 
itly in the reference. 

Decomposition of a Majoranon into two Majo- 
rana fields - Here, it is useful to go back to Eq. ([5]) 
briefly and discuss an interesting fact. Using the uni- 
tary operation U = ze"*'^'^"/'' ® ^-i-^iy^/i -^^ith the basis 
vEf = f7$ e g^p]^ ^Yig^^ ^ ^ MU^, the "Majorana 
Hamiltonian" , Hm in Eq- (0) becomes 

H = U^HmU, (13) 

= (I2 ® + "T-(ct2 (g) Cr^). (14) 

This Hamiltonian leads to two uncoupled Dirac equations 

idt^± = H±(l)±, (15) 
H± = axPx±rnaz, (16) 

with <j)+ = ($(i),$(2)) and 0_ = ($(3)^ $(4)) g^j,j^ ^Yiat 
^ = ((/)+,(/)_). It is interesting that are invariant 
under charge conjugation, which in 1-t-l D is defined as 
—icTzCyc/jj-. The uncoupling thus suggests that a Majo- 
ranon in 1-1-1 D can be decomposed into two Majorana 
fermions. In fact it is well-known that a field (which in 



1- 1-1 D has two components) obeying the Dirac-Majorana 
equation, a general relativistic field equation that has 
both the Dirac and the Majorana mass terms, can be 
broken down into two majorana fermions with different 
masses [l^iliBl- We provide a simple proof here. 

To construct a general proof for the Dirac-Majorana 
equation, we first write it in the Majorana representa- 
tion (in which charge conjugation is equal to complex 
conjugation): 

t'j'^df^'^ = mA/** -f- mn^'. (17) 

Decomposing the four component spinor as = (4*+ + 
i5'_)/-\/2 where are the Majorana fields, i.e., = 
it is readily seen that 

Z7^a^*± = imn ± rriM) *±. (18) 

Therefore a Dirac-Majoranon obeying Eq. ([TT]) can be 
decomposed into two Majorana fermions obeying their 
respective Dirac equations. The following proposal for 
quantum simulation thus work not only for a Majoranon 
but a Dirac-Majoranon too. 

A proposal for quantum simulation - The decom- 
position described above has an interesting consequence 
for quantum simulation of a Majoranon: it is possible to 
simulate the Majorana equation by simulating two Dirac 
equations with the opposite mass signs. An advantage is 
that only a single qubit is required, so the single-qubit 
addressibility and qubit-qubit interaction is not an issue 
anymore. The price to pay is that the method requires 
complete state reconstruction, including the continuous 
degree of freem. This is because the full information 
on spinors, (f)^ and is required in order to calculate 
expectation values. Despite the caveat, the method pro- 
vides a good alternative as far as observing exotic physics 
goes, especially for the case that only requires qubit state 
reconstruction, which we describe below. Note that the 
exotic physics of the Majorana equation does not man- 
ifest itself in the relativistic regime, in which the mass 
term becomes negligible and the Majorana equation re- 
duces to the Dirac equation. Instead it is most prominent 
in the stationary case as discussed in [l3l ■ 

It is also possible to simulate the TCME in 2 spatial 
dimensions, i.e., simulate a two-component chiral field in 

2- 1-1 D, using the fact that a Majoranon can be decom- 
posed into two Majorana fermions. To see this, first note 
that the (right-chiral) two-component Majorana equation 
and the Majorana equation are equal and can be written 
as 

idtip = {cTxPx + cryPy)ip " imay^}* . (19) 

By decomposing 1/; = ('!/;_|_-t-i'i/;_)/\/2 and assuming ili± = 
—iazCTy'iJj^, it is easily seen that the two Dirac equations 

idti>± = -{o'xPx + cryPy)'ip± ± m(Tzip±, (20) 

are equivalent to Eq. (jl9p . The same decomposition does 
not work in 3-1-1 D because the charge conjugation op- 
erator cannot be constructed for two-component spinors. 
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This seems to corroborate the fact that chiral spinors 
are the irreducible representations of the Lorentz group, 
i.e., they cannot be decomposed into smaller representa- 
tions. One is forced to use the four-component spinors in 
this case, requiring a quantum simulator in 3+1 D pro- 
posed in [3l for example. An interesting open question is 
whether two qubits instead of a single four-level system 
can be used to simulate the Dirac equation in 3-1-1 D. 

Simulation procedure - Thus far, we have shown 
that in 1-1-1 or 2+1 D a (two-component) chiral spinor, 
obeying the TCME or equivalently the Majorana equa- 
tion, can be either turned into a four-component Majo- 
rana field or decomposed into two two-component Majo- 
rana fields. The latter allows quantum simulation of the 
Majorana equation in 1+1 (2+1) D with a single qubit 
and one (two) continuous degrees of freedom, given that 
one can reconstruct the states completely. The simula- 
tion procedure can be summarized into four steps: (z) 
Choose a system that simulates the Dirac equation; e.g., 
trapped ions [1, 3, stationary light polaritons [l^ or op- 
tical lattice [i^ schemes, (ii) Prepare the initial state 
that you want to simulate; i.e. write down the spinor 
corresponding to a particular initial state and decom- 
pose it into two spinors obeying the Majorana conditions: 
if,' = + i'ip-, where (V'±)c = V'i- We will show specific 
examples below, [iii) Time evolve the states il)± accord- 
ing to the Dirac equation with ±to, respectively, (iv) 
Evaluate an observable in terms of ■0+ and tp- ■ 

In 1+1 or 2+1 D, the Majorana condition x± = 
—iazCryX± is satisfied by the basis vectors x+ ~ 
(1,-1)/V2 and x- = {i,i)/V2, and one can expand any 
initial state in terms of these basis spinors. A given state 
-0(0) can be expanded as {^^{0) + i'ilj-{0))/V2, in terms 
of the Majorana fields i/'i using the following relations 

V'+(O) = i[V;(O)-ia,a,0*(O)], (21) 

For example, ip{0) ~ (1,0) is equivalent to ip+{0) = 
(1,-1)/V2and V'-(O) = -{i,i)/V2. This example corre- 
sponds to a particle at rest discussed in [T^] . In this case, 
there is no momentum dependence and the simulation 
only requires qubit state reconstruction j20|. Therefore, 



one can simulate the Majorana equation in the rest-frame 
with a single trapped ion instead of two. As explained 
above, this corresponds to the limiting case where the dif- 
ferences between the Majorana and the Dirac dynamics 
are most prominent. 

As another example, consider a general state that dis- 
tinguishes between the Majorana and the Dirac equation, 
i.e., a state that is not charge conjugate invariant, 

^(0) = e'fo^e"^ Q . (23) 

This state can be created in the trapped ion setup . It 
is decomposed into 

■0+(O) = \/2 sin(poa;)e~4A2;(;_, (24) 

ip-{0) = ~V2 cos{pox)e~4A^x-- (25) 

In this case, evaluation of an observable requires com- 
plete information on 'i/'±j including the spatial depen- 
dence. This means the full quantum state tomography of 
the qubit plus the phonon mode for trapped ions, whereas 
in other setups, such as stationary light polaritons or an 
optical lattice, it would mean a spatially resolved detec- 
tion scheme. 

Conclusion - A Majoranon, i.e., a particle that obeys 
the Majorana equation but not the Majorana condition, 
can be decomposed into two Majorana fermions, i.e., par- 
ticles that answer to both the Majorana equation and 
condition. It is therefore possible to simulate a chiral 
spinor with previously proposed quantum simulators of 
the Dirac equation. Quantum simulation in 1+1 or 2+1 
dimensions requires a single qubit plus one or two con- 
tinuous degrees of freedom, instead of two qubits as pro- 
posed before, with the trade-off that complete state re- 
construction is required. An interesting exception is for 
a particle initially at rest, for which the differences be- 
tween the dynamics of the Majorana and the Dirac equa- 
tions are most obvious. In this latter case, only single 
qubit state reconstruction is needed, greatly simplifying 
the complexity of a possible experimental implementa- 
tion of our approach. 
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